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Abstract 

Energy loss of a jet due to multiple collisions in the nuclear target is studied in the framework of the 
perturbative QCD with an infrared cutoff. An iterational procedure for its calculation is developed, 
• which allows to reliably find up to 16 successive collisions in rescattering. The calculated shift of 

the scaling variable reaches values of the order 0.1 at medium x and transverse momentum p in the 
interval 5 — 10 GeV/c. The shift is found to be independent of the atomic number of the target and 
^ ■ energy. 
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\^ , Jet quenching in nuclear matter due to gluon emission has long been a subject of detailed studies, 

■ related to the observed behavior of particle spectra in heavy-ion collisions at RHIC [1] . Much attention 

has also been devoted to a simpler quenching mechanism due to hard elastic collisions ("collisional 
quenching"), which is formally of the leading order in as compared to gluon emission corresponding 
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' to inelastic hard collisions. It is expected that quenching due to elastic collisions is proportional to the 

■ length L travelled by the jet inside the nucleus, whereas gluon emission seems to lead to quenching 

proportional to [2]. However this conclusion has been derived for superhigh energies. The relative 
strength the two sources of quenching at present energies is not clear both from the theoretical and 
experimental points of view [3]. 

Most calculations of the collisional energy loss (GEL) refer to the propagation of a fast parton 
in the homogeneous quark -gluon plasma [4, 5, 6, 7]. The study of parton energy loss in a realistic 
nucleus is a harder task. It requires using the multiple scattering formalism beyond the leading 
Glauber approximation with full attention to the interplay between the longitudinal and transverse 
momenta. First calculations made in [8] indicated that elastic collisions diminish the jet energy quite 
considerably, up to 40% at LHG energies. However due to technical difficulties the authors of [8] 
were able to push their calculations of multiple hard scattering of the jet inside the nucleus only up to 
the first three collisions. This is quite insufficient even for RHIG energies at which, as our calculation 
show, at least 6 collisions are to be taken into account. At the LHG energies this number grows to 
12. To obtain reasonable results for the overall quenching due to elastic collisions one has to search 
for methods which allow to study multiple parton collisions with the loss of energy with a reasonable 
accuracy. This paper presents an iterative approach which seems to give satisfactory results at both 
the RHIG and RHIG energies. 
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Our framework follows that of [9, 10, 11] which considers multiple scattering of a parton existing 
inside the fast projectile long before the actual collision. The scattering itself is due to the QCD 
perturbative interaction in the lowest order, with a relatively high transferred transverse momentum 
and a loss of energy determined by the kinematics. Naturally we introduce an infrared cutoff at low 
transferred momenta (by an effective gluon mass) and our results depend on the value of this cutoff. 
On physical grounds we choose it to lie in the region of 1.5-^ 2.0 GeV/c. With such values our iterative 
procedure gives reliable results up to 16 collisions at the LHC energy. 

From the start we have to stress that we consider only passage of a jet through the nucleus. The 
study of its hadronization is postponed for future publications. We study the cross-sections as a 
function of the transverse momentum p of the observed jet and its scaling variable x. For the study 
of collisional energy loss of a jet due to rescattcring in the nucleus use of these variables is natural. 
Hadronization, necessary to study the observed parton spectra, can be conveniently calculated in these 
variables. To pass to c.m.rapidities instead of x one can use the standard relation y = ln{xW/p) where 
W is the cm. energy for proton-proton collisions. 

The CEL effect was found to be quite noticeable. For values of p in the interval 5 — 10 GcV/c 
and X of order 0.5 the decrease in the scaling variable was calculated to lie in the interval 0.09 — 1.14. 
However a curious result is that this decrease does not depend on the energy nor on the atomic number 
of the target. The latter implies that CEL does not depend on the length of the rescattering path 
inside the nucleus. Cross-sections themselves exhibit more or less the expected behaviour: they fall 
with the inclusion of CEL. This fall is only weakly energy dependent, so that the general picture at 
the RHIC and LHC energies remains practically the same. As a curiosity at the RHIC energy we 
found that the standard nuclear suppression factor Ra at small values of x actually grows with the 
introduction of CEL. This abnormal behaviour however is not seen at higher values of x nor at the 
LHC energies. 

2 Glauber-like parton rescattering. 

To formulate our model and have a benchmark for the following calculations we briefly reproduce 
here the derivation of the inclusive jet production in hadron-nucleus or nucleus-nucleus scattering in 
the pure Glauber approximation without loss of energy. The details may be found in the original 
derivation in [11]. 

We assume the nuclear state |A > to be a superposition of states with a different number n of 
partons, each characterized by its scaling variable x and impact parameter b combined into argument 
z = {x, b}: 



n i=l 

If the observed parton originates from the A-nucleus (moving along the z-axis) , then the wave function 
of the latter is for a given n: 



Here V'a(-2^) denotes the wave function of the observed parton, and a = {x,p}, where x and p - are the 
scaling variable and transverse momentum respectively. 




(1) 



Zn)- 
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The inclusive cross-section of parton production at a given impact parameter (3 for the nucleus- 
nucleus (AB) collision is given by the following expression 

do^ f 

^ = Z]"' / dzidz[i^a{z[)i^a{zi)dTA{n-l)dTB{l)'^Xn-li^l^Z2,-,Zn) 
^ nl 

x'^A,n-i{zi,:,Zn)\'^B,i{uj)\^[S*iiz'i,:,Zn\ui,:,ui) - l][Sniizi,:,Zn\ui,:,ui) - 1], (3) 

where dT{n) is the pase volume for the state with n partons. This expression contains contributions of 
both hard and soft parton interactions. The latter interactions are disregarded because they cannot be 
studied in the perturbative QCD. The contribution of hard collisions is given by the product S*; x Sni- 
Further transformations rely on the assumption of factorization of nuclear S-matrix into a product 
of elementary S-matrices corresponding to parton-parton interactions, which is the essence of the 
Glauber model: 

n / 

Snl{zi,...,Zn I Ui,...,Ui) = 11 W 

i=lj=l 

Here Sij is S-matrix for the interaction of parton i of the A nucleus with parton j of the B nucleus, 
while Uj - refers to parton variables of nucleus B. We assume factorization of parton distributions: 

l*A(B),Kn,-)|^ = ^e-wnrK), (5) 

with (/) = / d?uTB{u) and Vb{u) denoting the product of the parton distribution inside the nucleon 
and nuclear profile function, 

^A{B) {u) = Ta{b) {c)Pa[B) {^)- (6) 

As a result, we get the final expression for the inclusive cross-section of jet production in AB collisions 
with hard rescattering taken into account: 

daAB 



where 



(27r)2_^i^ = / dHd\e'^TA{b - ^)P^(x)e-^«(^)^«(-'°){e^« - 1}, (7) 

d^pdxd^p J 

FB{x,b) = J ^ J dujPB{u)I{x,uj,q)e'P' (8) 

and I{x,LO,q) is the distribution in the transverse momentum q in the elastic collision of two partons 

with their scaling variables x and lo. 

For hadron-nucleus collisions one has to change T into (5-function and integrate over /? to obtain 

(2vr)2^i = I d2MVe^^"-p(x)e-^^(^)^^(-''^){e^-*W^^(-''-) - l}, (9) 

where P(x) is the parton distribution in the incoming hadron. 
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3 Energy loss in a single elastic parton collision 

3.1 Kinematics 

Following the picture of multiple parton scattering let an incoming (massless) parton (4-momentum 
p = {p+,p-,p±)) scatter on a target parton (4-momentum k = {0,k-,0)). The final momenta are 
p' = {p'_^_,p'_,p± — q) and k' = {k'_^, k'_, q) where q is the transverse part of the transferred momentum. 
The scaling variables of the initial and final projectile parton are x and x', those of the target parton 
w and w' . The parton collision cm. energy squared is sxw where s is the overall cm. energy squared. 
The longitudinal momenta conservation gives 

The energy loss of the incoming jet is determined by the difference 

Ap+ =P+- P+ = (11) 

which transforms into the corresponding decrease of the jet scaling variable x. 
Prom Eqs. (1) we conclude 

2k'_p'+ + q^ = 2k'_p+, 2k'_p'+ + {p±-qf = 2p'^^^ + k.^ (12) 

and thus 

q' - 2kLp+ = (px - qf - 2p'+ + k) . (13) 



This allows to express k'_ via p'_^ : 



k'_ 
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2p+l'^'2p+ 

and obtain an equation: 



2p'i^ + k.) + q^-ip^-qf] (14) 



(^^ - 1 j (2pV(^ + k.) +q'- {p^ - qf) +q' = 0. (15) 

This is a quadratic equation from which one determines the initial projectile longitudinal momentum 
p^ in terms of its final longitudinal momentum and other variables describing the collision: w, q and 
p±. Denoting 

one finds from (15 

u = + P±q - V iP+k- + P±q)^ - pW}- (17) 

P± 

For small transferred momenta u is proportional to q^ and then 

q'^ 
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so that ^ 

^ = 277 A (19) 

1 — q'^/[swx') 

and is independent of p±, although the latter may be comparatively large as a result of accumulation 
of several transferred momenta. We use this simplified formula in our calculations. 

Note that condition < x < 1 leads to the restriction q"^ < sw{x' — x'^), which limits the region of 
integration in g^, x' and w. Prom this inequality it also follows that in any case < s{x' — x''^) < s/4. 

3.2 Parton-parton cross-section 

In principle one should take into account partons of different flavor both in the projectile hadron 
and target nucleus. Although it is straightforward, it considerably enhances the calculation time. To 
economize on the latter, following [12], we considered the effective gluon jet defined by the density 

P{x, Q2) = g(a:, Q2) + 1 (^q(x, Q^) + q{x, Q^)) , (20) 

where g,q and q are the gluon, quark and antiquark densities respectively. Correspondingly for the 
parton interaction we took the gluon-gluon one with the cross-section, integrated over the target gluon 
distribution [12] 

'^^ = dx'd^q = ^^""^ W J ^[^^(^)]-'+<^V-'[^G'H]^,+,2/,,(l - , (21) 

Here, according to our notation, x and w refer to initial gluons and x' and w' refer to final gluons. We 
pass to the integration over the initial scaling variable of the target 

I{x',q) = (27r)2|4 / ^^[xG{x)],,+^y,^,wG{w)(l - -^f . (22) 

2q^ J w — q-^/sx' ^ sxw^ 

Prom (16) and (18) one finds 

w — q^ / sx' = wil — u), 1 = 1 — u + u^, 

sxw 

so that the cross-section becomes 

,9a^ f dw (1-a + o^)^ 



I{x',q) = (27r)2|^ J ^xG{x)wG{w) 



2q^ J w ^ ' ' ' I- a 

~ (27r)^^ J ^xG{x)wG{w){l - uf (23) 

where the second approximate form takes into account that q^ and consequently u arc assumed to be 
small. The cross section I{x' , w, q) entering our formulas for the n-fold cross-section is obtained from 
this one by dropping the partonic distributions and integration over w. 
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4 Loss of energy 



Eqs. (7) or (9) have an obvious physical interpretation. Developing the exponents in powers of F one 
gets, say, for hA collisions the expression for the n-fold rescattering 

T (n) 1 /• 

(27r)2^^ = - J dHd^re'P^P{x)T2{b){[FA{x,r) - FAiO,r)y - F:i{0,r)} (24) 

At fixed r one finds the product of n inclusive cross-sections FA_{x,r) summed with similar products 
in which some of the inclusive cross-sections are substituted by the total cross-sections Fa{x,0) with 
a minus sign. The second term in (24) eliminates the product which does not contain inclusive cross- 
sections at all. In the momentum space the product goes over into a convolution. So physically this 
term corresponds to a number < A; < n of consecutive inclusive cross-sections and (n — k) total 
cross-sections distributed among these inclusive cross-sections in all possible ways. 

This indicates how one can take into account the loss of energy during rescattering. Each inclusive 
cross-section in our approximation corresponds to the elastic scattering with a change of transverse 
momentum. If after n such collisions the observed scaling variable of the jet is x then before the 
n-th elastic collision and so after the (n — l)-th collision the jet has its scaling variable Xn-i > x. 
Similarly before the (n — l)-th collision and after the (n — 2)-th collision the jet has its scaling variable 
Xn-2 > Xn-i and so on. In the first elastic collision the scaling variables before and after the collision 
are xq and xi respectively. If the loss of energy would be independent of the scaling variable of the 
target parton w then this loss would be accounted for by a simple substitution 

n 

F2{x,r)^l[F{xi,r) (25) 

i=l 

where Xi arc the scaling variables of the jet after the i-th clastic collision with Xn = x the scaling 
variables of the observed jet. However, as we have seen, the change in x in fact depends on w, so that 
this substitution has to be done under the sign of the integration over all scaling variables Wi of the 
target partons. 

Another problem is to distribute the total cross-sections among the clastic ones with correct scaling 
variables. Obviously the total cross-sections appearing between the A;-th and (k + l)-th elasic ones 
are to depend on Xk- For instance for n = 2, with the loss of energy independent of w, the original 
Glauber expression 

F^(x, r) - 2F{x, r)F{0, r) (26) 

should be substituted by 

F{xi,r)F{x2,r) - F{x2,r)F{xi,0) ~ F{x2,r)F{x2,0), X2 = x (27) 

However the -w-dcpcndcnce again requires making this substitution inside the integral over all w's. 

A convenient way to obtain constructive formulas which take into account the loss of energy is to 
set up an iteration procedure allowing to find the (n -|- l)-fold rescattering contribution in terms of 
the n-fold one. 

4.1 Iteration procedure 

It is more convenient to work directly in the momentum space. Consider a sequence of n collisions, 
some elastic and some total. Call the contribution of the amplitude which ends with the n-th elastic 
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collision, with nuclear factor separated, X„(a;,p) where x andp are the scaling variable and transverse 
momentum of the jet. Then the total contribution with any number of elastic collisions after the last 
inelastic one is 

n 

Y^"\x,p) = Xn-k{x,p)[-a{x)]\ (28) 

fc=0 

where cr(x) is the total cross-section of the projectile parton and by definition 

Xq{x,p) = {27rf 6\p). (29) 

Prom (28) a relation follows 

yr(^,0 = Xnix,r) - C7(a;)y„*°*i(a;,r). (30) 

Note that Y^°^{x,r) includes an improper term corresponding to the case when there are only elastic 
collisions: 

rr(x,p) = Y^{x,p) + [-a{x)r{27Tf6\p). (31) 
Here 1^ is a contribution from n collisions with at least one inelastic. One easily finds that also 

Yn{x,p) = Xn{x,r) - a{x)Yn-i{x,p). (32) 

Our basic equation connects Xn with Under the sign of integration over the target scaling 

variable w we have just to convolute Y^f'{xn-i,p) with F{xn,p) to obtain 

Xn{x,p) = J ^dwnP{w)I{q,x,w)Y^°MxuP-q)- (33) 

Here x is the scaling variable of the projectile AFTER the n-th collision, q and w are the transferred 
momentum and the scaling variable of the target in this last inelastic collision. The scaling variable 
xi is the one BEFORE the n-th inelastic collision considered as a function of collision variables 

xi=f{q,x,w). (34) 

It is generally greater than x, which corresponds to collisional quenching. The explicit expression of 
function / is given by (19). The momentum p — q is the projectile momentum BEFORE the last 
inelastic collision. With only one elastic collision 

Xi{x,p) =) = y ^dwP{w)I{q, X, w){2'K)H\q - p), (35) 

which corresponds to 

Y^"\x,p) = {2nfS\p) (36) 

and follows from our assumption that the projectile has no transverse momentum before all collisions. 
Separating the improper contribution to in (33) we get 

Xn{x,p) = j ^2^cZt(;P(«;)/(g,x,u;){y„_i(xi,p-g) + (-(7(xi))""'(27r)2,52(p-g)} (37) 
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and using relation (32) get the final form of the iteration equation 



Yn{x,p) = J ■^^dwP{w)I{q,x,w)^^Yn-i{xi,p - q)+ 

( - a{xi)y~\27r)^6\p - q)} - a(x)y„_i(x,p) 
One starts iterations with 

Yo{x,p) = 

and successively determines Ti, T2 and so on. At fixed b the final contribution to the inclusive 
section is obtained as 

This iterative procedure can be directly rewritten in the r space. One obtains 

/(Pq 
j—^e"'''dwP{w)I{q, x, w)Yn-i (xi {q, x, w), r) 

/ci^q ■ ( \ra— 1 

-^—^e'-'i"' dwP{w)l{q,x,w)\^ - a{xx{q,x,w)) j - a{x)Yn-\{x,r). 

4.2 The master equation 

Introduce 

z„ = c\TmYr,{x,p) ~ l(^r^(6))-y,(x,p) = 

where we denoted ^ = ATa^^)- Obviously 

-,e-^Yn{x,p). 



dZn 1 ^n-U 



(n-iy. 

Now take the recurrent relation (38) and multiply it by /{n — 1)! to get 
dZn f d^q 



—--^dwP{w)I{q,x,w)^Zn-iixi,p - q) 



J (27r) 

+-^^^^( - <j{xi)y~\27r)H\p - q)} - a{x)Zn-i{x,p). 
Summing from n = 1 to 00 and defining 

00 

Z{C,X,p) = Zn, 
n=l 

SO that 

^ ' - P{x)Z{^,x,p), 



dxd^pd^b 
we obtain a master equation for Z 



+e-^''^'^'\27Tf6\p - q)} - a{x)Z{C,x,p). (47) 

Technically this equation does not present serious advantage as compared to the direct iteration 
procedure, since it allows to calculate the contribution only at fixed impact parameter b and requires 
subsequent integration over b, which for a fixed number of collisions can be done more easily. 

5 Numerical results 

As mentioned, the iteration procedure allows to reliably calculate rescattering contributions up to the 
16th successive collision. We choose the gluon mass m = 2 GeV/c as an effective infrared cutoff. 
So our inclusive cross-sections are meaningful at p± > 2 GeV/c. We study hA collisions with gold 
as a target {A = 197). For comparison also the copper target was studied (A = 64). The standard 
parameter K, which effectively takes into account higher orders in parton scattering cross-sections has 
been chosen to be 1.5 at the RHIC energies and 1.0 at the LHC energies. For the parton distributions 
we used the CTEQ6L set [13]. 

At the RHIC energy {^/s = 200 GeV) the rescattering cross-sections fall very rapidly with the 
number n of the subsequent collisions and are practically zero at n > 6. In contrast at the LHC 
energy (-/i = 8800 GeV) the decrease of the cross-sections with n is much slower and even the 12th 
collision contributes 2 5% of the total cross-section. 

Our results for the RHIC energy are presented in Figs. 1-6. In Fig. 1 we illustrate the effect of 
collisional energy loss (GEL) presenting the ratio of the cross-section with and without GEL 

Q{x,p) = ^/(^) (48) 

dxd'^pi ^dxd'^p' no energy loss 

(In the cross-section with no energy loss x is conserved during rescattering and is equal to x of the 
observed jet). As we observe the rate of CEL grows both with x and p, somewhat flattening at p > 10 
GcV/c. In Fig. 2 we compare absolute values of the cross-sections with and without CEL as a function 
of X at two values of p = 5 and 10 GeV/c of the transverse momentum. These curves serve us to 
crudely determine the magnitude of the energy loss A as a difference between the values of x at which 
cross-sections with and without energy loss give the same cross-section Thus found values of A are 
presented in Fig. 3 as a function of x ^. They show a maximum at medium values of a; ~ 0.5 with 
A = 0.09 for p = 5 and 0.14 for p = 10. The growth of collision loss with p was to be expected from 
Eq. (19). 

An unexpected feature of the CEL is its practical independence of the atomic number of the target. 
In Fig. 4 we compare the shift A of the scaling variable for Au and Cu at p = 5 GeV/c. The two 
curves practically coincide. In the common language it implies that the energy loss does not depend 
on the length of the path inside the nucleus at all. This can be understood by the fact that the optimal 
configuration of rescattering to give the observed transverse momentum is to produce this momentum 
in just one collision, leaving the rest collisions to have a small momentum transfer. Then obviously it 
does not matter how many scattering centers there are along the rescattering path. 

Finally in Figs. 5 and 6 we show the nuclear suppression factor (NSF) 



^The chisen primitive procedure to find A fails t high values of x where the cross-sections become very small, which 
explains irregularities of the curves at the right end. 
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Figure 1: The ratio of the cross-sections with colhsional energy loss to the ones without it (Eq. 48) for 
p-Au colHsions at -y/i = 200 GeV. The curves from top to bottom correspond to x = 0.025, 0.2, 0.4, 0.6 
and 0.8. 

as a function of p two values of x = 0.025 and 0.4 with and without CEL. The curious result is at low 
values of x introduction of CEL actually enhances the NSF ('antiquenching effect'). This is partially 
the result of the above mentioned optimal configuration for rescattering in which the effect of energy 
loss may be stronger for the single rescattering contribution than for the total one. 

Our results for the LHC energy ^/s = 8800 GeV and Au target are presented in Figs. 7- 10. 
Fig. 7 illustrates the CEL showing ratio Q{x,p) (Eq. (48)) as a function of p at different values 
of X. The general picture remains the same as at the RHIC energy (cf. Fig. 1). However the x 
dependence becomes stronger, so that quenching at low x becomes weaker and at larger x the same 
or even stronger. In Fig. 8 we compare the shift A of the scaling variable at p = 5 GeV/c at LHC 
and RHIC energies. There is practically no change at all, so that according to our calculations CEL 
not only is A independent but also energy independent. Finally Figs. 9 and 10 show the NSF as a 
function of p at x = 0.025 and 0.4 respectively. At x = 0.025 the effect of CEL is much weaker than 
at the RHIC energy and acquires the expected direction: it lowers the NSF. At re = 0.4 CEL lowers 
the NSF much stronger than at the RHIC energy. 

6 Conclusions 

We have calculated colhsional energy loss of emitted jets due to multiple collisions of partons from the 
projectile hadron with partons in the target nuclei. Our results indicate that the loss grows with the 
transverse momentum p of the observed jet and is maximal at medium values of the scaling variable 
X, where the shift in x reaches values in the interval 0.09 — 0.14 for p varying from 5 to 10 GcV/c. A 
somewhat unexpected result is that this shift is practically independent on the atomic number of the 
target nor on the energy, being the same for RHIC and LHC. The ^-independence implies that the 
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Figure 2: Cross-sections for p-Au collisions at ^/s = 200 GeV with and without collisional energy loss 
as a function of x for p = 5 GeV/c (upper pair of curves) and p = 1^ GeV/c (lower pair). In each pair 
the upper curve corresponds to the absence of energy loss. 
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Figure 3: The shift in the scaling variable at p = 5 GeV/c (lower curve) and 10 GeV/c for p-Au 
collisions at = 200 GeV. 
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Figure 4: The shift in the scahng variable at p = 5 GeV/c for p-Au (the sohd curve) and p-Cu (the 
dashed curve) colhsions at ^/s = 200 GeV. 




Figure 5: Nuclear suppression factor for p-Au collisions at -y/s = 200 GeV as a function of transverse 
momentum p at x = 0.025 with (the solid curve) and without (the dashed curve) collisional energy 
loss. 
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Figure 6: Same as Fig 5 at x = 0.4. 
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Figure 7: The ratio of the cross-sections with colhsional energy loss to the ones without it (Eq. 48) for 
p-Au coUisions at ^/s = 8800 GeV. The curves from top to bottom correspond io x = 0.025, 0.2, 0.4, 0.6 
and 0.8. 
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Figure 9: Nuclear suppression factor for p-Au collisions at -y/s = 8800 GeV as a function of transverse 
momentum p at x = 0.025 with (the solid curve) and without (the dashed curve) collisional energy 
loss 
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Figure 10: Same as Fig 9 at x = 0.4. 

shift docs not depend on the length passed by the jet inside the nucleus. An explanation may be that 
most of the shift takes place in just one collision during rescattcring. 

We have also calculated the so-called nuclear suppression factor (NSF) with collisional energy 
loss taken into account. From our results it follows that at moderate transverse momenta NSF still 
stays noticeably above unity, indicating that energy loss does not reverse the enhancement in the 
j5-distribution due to kicks in collisions with partons in the target. Surprisingly at the RHIC energy 
and small x the energy loss even enhances NSF, which again may be the consequence that the loss is 
mostly due to just one collision. However one should take into account that our treatment considers 
the loss in a single collision to be relatively small, which condition may be violated at small values of 
xs (see Eq. (19)). 

Finally we stress that our results refer to production of jets. To pass to observed hadrons one 
should convolute our cross-sections with the appropriate fragmentation functions. Although this 
procedure seems direct, we prefer to postpone it for future studies, because it inevitably involves new 
badly controlled assumptions and approximations, especially when hadronization takes place inside 
the nucleus. We prefer to present results for jets, strictly within the perturbational QCD framework, 
which do not depend on these soft physics assumptions. They may constitute a starting point for the 
study of hadronization process. 
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